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Abstract: The partition function on the three-sphere of many supersymmetric Chern-Simons- 
matter theories reduces, by locahzation, to a matrix model. We develop a new method to study 
these models in the M-theory limit, but at all orders in the expansion. The method is based 
on reformulating the matrix model as the partition function of an ideal Fermi gas with a non- 
trivial, one-particle quantum Hamiltonian. This new approach leads to a completely elementary 
derivation of the N'^l'^ behavior for ABJM theory and = 3 quiver Chern-Simons-matter 
theories. In addition, the full series of 1/N corrections to the original matrix integral can be 
simply determined by a next-to-leading calculation in the WKB or semiclassical expansion of the 
quantum gas, and we show that, for several quiver Chern-Simons-matter theories, it is given 
by an Airy function. This generalizes a recent result of Fuji, Hirano and Moriyama for ABJM 
theory. It turns out that the semiclassical expansion of the Fermi gas corresponds to a strong 
coupling expansion in type IIA theory, and it is dual to the genus expansion. This allows us to 
calculate explicitly non-perturbative effects due to D2-brane instantons in the AdS background. 



Contents 



1. Introduction 1 

2. The ABJM matrix model in the 't Hooft expansion 5 

2.1 1/N expansion and non-perturbative effects 5 

2.2 The partition function as an Airy function 6 

3. Chern— Simons— matter theories as Fermi gases 9 

3.1 ABJM theory as a Fermi gas 9 

3.2 More general Chern-Simons-matter theories 14 

4. Thermodynamic limit 17 

4.1 The thermodynamic limit of ideal Fermi gases 17 

4.2 A simple derivation of the N^/"^ behavior in ABJM theory 18 

4.3 Large N corrections 20 

4.4 Relation to previous results 22 

5. Quantum corrections 22 

5.1 Quantum-corrected Hamiltonian and Wigner-Kirkwood expansion 22 

5.2 General structure of quantum corrections 25 

5.3 Quantum corrections in ABJM theory 27 

5.4 Quantum-mechanical instantons as worldsheet instantons 32 

6. More general Chern Simons matter theories 34 

6.1 Thermodynamic limit for general necklace quivers 35 

6.2 Airy function behavior for a class of CSM theories 37 

6.3 ABJM theory with fundamental matter 40 

6.4 The massive theory 43 

7. Conclusions and prospects for future work 46 
A. Quantum corrections in ABJM theory at order 0{ti^) 48 



1. Introduction 

One of the most interesting aspects of the AdS/CFT correspondence is that, in principle, one can 
use gauge theories to learn about elusive aspects of string theory and quantum gravity. For ex- 
ample, ABJM theory [1], as well as other supersymmetric Chern-Simons-matter (CSM) theories, 
are conjecturally dual to M-theory on spaces of the type AdS4 x Xy. Therefore, computations in 
the gauge theory side might give interesting insights on M-theory on these backgrounds. 

It was shown in [2] that the partition function on the three-sphere of CSM theories with 
> 3 supersymmetry reduces, via localization, to a matrix model. This result was extended to 
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theories with M > 2 supersymmetry in [3, 4]. In the case of ABJM theory, the corresponding 
matrix model was solved for arbitrary 't Hooft coupling and at all orders in in [5], providing 
the first gauge theory derivation of the famous N^/"^ growth of the number of degrees of freedom 
of N M2 branes [6] . In the last year, many results have been obtained for these matrix models, 
providing precision tests of the AdS4/CFT3 correspondence as well as beautiful field-theoretical 
results on supersymmetric CSM theories (see [7] for a review and a list of references). 

Up to now, the "stringy" side of these matrix models has been studied less intensively, but 
there have been already some interesting results in this direction for the ABJM theory. In [5] 
it was found that the genus g free energies of the matrix model contain very rich information 
about worldsheet instantons in type IIA theory (i.e. they include non-perturbative effects in a'). 
In [8], by studying the large order behavior of the genus expansion, it was possible to identify 
as well non-perturbative corrections in the string coupling constant, conjecturally associated to 
D2-branes, or to membrane instantons in M-theory. Finally, building on the results of [5, 8], it 
was shown in [9] that, once worldsheet instantons are discarded, the full free energy of the ABJM 
matrix model is given by an Airy function. Schematically, we have 

ZAB3M(xAi[f{k){N-g{k))], (1.1) 

where k is the CS level (i.e. the inverse coupling), f{k) oc k^^^ is a function of k determined by the 
large N limit, and g{k) is a function of k which shifts N (see (2.14) below for a precise formula). 
This is a beautiful result which provides the all-orders expansion of the partition function in 
powers of the string length, and therefore resums the perturbative long-distance expansion for 
quantum superstrings. 

This result raises an interesting possibility. The N^/"^ behavior of the free energy characterizes 
a large class of A/" = 3 [10, 11] and J\f = 2 [12, 13, 14, 15] Chern-Simons-matter theories. On the 
other hand, the result (1.1) says that, for ABJM theory, this behavior is just the leading term in 
the logarithm of the Airy function. It is then natural to propose the following 

Conjecture 1.1. In parity-invariant supersymmetric CSM theories that display an N^/"^ growth 
in the number of degrees of freedom, the leading large limit and the 1/A^ corrections to the 
partition function on the three-sphere add up to an Airy function, i.e. we have schematically 

ZcSM cx Ai [f{ka) {N - g{ka))] , (1.2) 

where ka are the different CS levels involved in the theory, f{ka) is a function which is determined 
by the leading, large N limit, and g{ka) is a shift which also depends on the details of the theory 
(but cannot be determined by the large N limit alone). 

If this conjecture is true, the corresponding matrix models display a universal behavior 
characterized by the Airy function. For theories which are not parity-invariant, we expect the 
Airy function to be the crucial ingredient of the answer. It is interesting to point out that the 
universal role of the Airy function in summing up 1/iV corrections was already proposed in [16], in 
the different but related context of five-dimensional black holes made out of M2 branes in certain 
Calabi-Yau compactifications^ . Unfortunately, the techniques used to derive (1.1) for ABJM 
theory rely heavily on the calculation of 1/N corrections based on the holomorphic anomaly 
equation (sec [5, 9] for details and references). For other models it is not clear how to generalize 
these techniques, even in cases (like the one studied in [11]) where the planar resolvent is known 
explicitly. 

^We would like to thank C. Vafa for discussions on this. 
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All the results mentioned so far have been obtained in what we will call the 't Hooft expansion 
of the matrix model, which corresponds to the genus expansion in type IIA superstring theory. 
The 't Hooft expansion is the asymptotic expansion as N goes to infinity and the 't Hooft 
parameter of the model, A, is kept fixed at large N, 

N 

N -^oo, A = — fixed. (1.3) 

However, in order to make contact with M-theory, we have to consider the M-theory expansion, 
i.e. the asymptotic expansion as goes to infinity in which k is fixed, 

N ^oo, k fixed. (1.4) 

The 't Hooft expansion of the matrix model gives some information about the M-theory expan- 
sion. For example, the all-orders result (1.1) presumably captures the all-orders expansion of 
the M-theory partition function in powers of the Planck length, at finite k. However, important 
contributions to the partition function in the M-theory expansion (like membrane instantons and 
Kaluza-Klein modes) are not directly captured in the 't Hooft expansion, and in order to use the 
ABJM matrix model as a tool to explore M-theory, one should study the regime (1.4) directly. 
A step in this direction was taken in [10], who found a very simple method to extract the large 
N, fixed k behavior of general CSM matrix models. However, in the approach of [10] it is not 
obvious how to compute systematically 1/N corrections to the leading large N behavior, not to 
speak about exponentially small corrections in N. 

It would then be very interesting to find a method to analyze the M-theory expansion of CSM 
theories, beyond the leading large N contribution considered in [10]. Such a method would allow 
us to address the above conjecture about the universal role of the Airy function in resumming 
the 1 /N corrections, and eventually could give us information about M-theoretic features of the 
matrix models which are not manifest in the 't Hooft expansion. 

In this paper, we make a first step in this direction, and we propose a new method to 
analyze the matrix model of some Chern-Simons-mattcr theories which fulfills some of the ex- 
pectations that we have just listed. The method consists of writing the partition function on 
the three-sphere as the partition function of an ideal Fermi gas with a non-trivial one-particle, 
quantum Hamiltonian. In this reformulation of the problem, the Chern-Simons level k becomes 
the Planck's constant h of the quantum-mechanical problem. Since k corresponds to the inverse 
string coupling, the semiclassical expansion of the Fermi gas is a strong coupling expansion in the 
type IIA theory large A'' dual. As usual, the large A^ limit is simply the thermodynamic limit of 
the gas. Our approach has the following features: 

1) The large A^ limit of the free energy at finite k is governed by the thermodynamic limit 
of the Fermi gas, which can be determined by a semiclassical calculation. In particular, we 
find a completely elementary derivation of the N^/'^ behavior of the free energy of ABJM theory, 
including the correct coefficient. This can be extended in a straightforward way to = 3 necklace 
quivers, and the result for the free energy is in full agreement with the calculation in [17, 18] 
based on the matrix model analysis of [10]. Interestingly, the relevant Fermi surface describing 
the thermodynamic limit of the Fermi gas is a two-dimensional polytope which characterizes the 
geometry of the dual tri-Sasakian spaces. In the case of ABJM theory, this Fermi surface is also 
a real version of the tropical curve obtained in [11]. 

2) The full 1/A^ expansion of the free energy of the matrix model is determined by the first 
quantum correction to the semiclassical limit. In this way we reproduce the result (1.1) for ABJM 
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theory at finite k, using again elementary methods in quantum Statistical Mechanics. Moreover, 
we prove our conjecture (1.1) for a large class of M > 3 supersymmetric CSM theories. 

3) One can also compute exponentially suppressed effects at large which are clearly M- 
theoretic. In particular, we find a systematic method to determine the contribution of membrane 
instantons. These non-perturbative effects receive however corrections at all orders in the h 
expansion, and so far we can only determine them order by order in k (but to all orders in 
the membrane winding). The strength of these corrections (i.e. the minimal membrane action) 
agrees with the instanton analysis of [8] . 

It follows from the last point above that, at the level of non-perturbative corrections, our 
method does not fully capture the M-theory expansion, since so far we are only able to determine 
these corrections in an expansion in k around k = 0. In order to make contact with the true 
M-theory expansion one should resum the resulting series. In spite of this limitation, the Fermi 
gas picture gives a concrete computational method to address non-pcrturbative effects in these 
superstring theory backgrounds. In fact, the semiclassical expansion of the Fermi gas is dual to 
the conventional genus expansion captured by the 't Hooft expansion of the matrix model. For 
example, in the 't Hooft expansion, worldsheet instantons appear as exponential corrections in 
the 't Hooft coupling, order by order in the Qg expansion, while membrane corrections appear 
as large N instantons, of order cxp( — l/(/s). In the Fermi gas approach developed in this paper, 
membrane instantons appear as exponential corrections in the chemical potential of the gas, 
order by order in the h expansion, while worldsheet instantons appear as quantum-mechanical 
instanton effects of order exp(— 1/^). 

Finally, we would like to point out that the use of Fermi gas techniques in the analysis of 
matrix models goes back to the solution of matrix quantum mechanics in [19], and should be 
familiar from the study of the c = 1 string (see for example [20] ) . The idea of studying the matrix 
integral partition function in the grand-canonical ensemble appeared in [21], and was developed 
in detail in [22, 23, 24]. In particular, the semiclassical limit of the Fermi gas was already used in 
Appendix A of [23] . The systematic application of semiclassical techniques of many-body physics 
in the study of these matrix integrals, which we develop in this paper, seems however to be new. 

The paper is organized as follows. In section 2 we review previous results on the 1/A'^ 
expansion of the ABJM matrix model in the 't Hooft expansion, focusing on [5, 8, 9]. In section 
3 we show that the matrix integral of a general class of M > 3 CSM theories (necklace quivers 
with fundamental matter) can be written as the partition function of an ideal Fermi gas with 
a non-trivial one-particle Hamiltonian. In sections 4 and 5 we present the tools to analyze the 
Fermi gas, and we illustrate them in ABJM theory. More precisely, in section 4 we study the 
Fermi gas in the thermodynamic limit, by passing to the grand canonical ensemble. This makes 
it possible to derive the leading N^^"^ behavior of the free energy of ABJM theory, by using 
elementary tools in Statistical Mechanics. We also compute exponentially suppressed corrections 
to the grand canonical potential, which are interpreted as membrane instantons in M-theory. In 
section 5 we study the quantum corrections to the grand canonical potential. We show that, 
up to non-perturbative terms, a next-to- leading WKB calculation is enough to determine the 
full 1/A^ expansion of the canonical free energy. This provides a simple derivation of the Airy 
function resummation of [9]. In section 6 we extend our techniques to more general CSM theories, 
including necklace quivers and theories with fundamental matter. We show that, when the free 
energy on the three-sphere is real, the 1/A^ expansion at fixed k gets resummed by an Airy 
function, thus proving conjecture 1.1 for this family of examples. We also consider the "massive" 
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theory of [53], where a different iV^/^ scahng has been found for the free energy, and we rederive 
it with our techniques. Finally, in section 8 we conclude with some prospects for future work. In 
an Appendix we collect some results for the grand canonical potential of ABJM theory at order 

As this paper was being prepared for submission, the paper [25] appear which also considers 
ABJM theory in the grand canonical ensemble. 

2. The ABJM matrix model in the 't Hooft expansion 
2.1 1/N expansion and non-perturbative effects 

The matrix integral describing the partition function of ABJM theory on is given by [2] 

■Z'abjm(-/V) 



1 



d^^ 



2 sinh (^^^) 



2 sinh 



2 / 



J (27r)^(27r) 



N 



n 



2 cosh 



exp 



\k_ 



N 



i=l 



(2.1) 



This matrix integral can be solved in the 't Hooft expansion (1.3) by using techniques of matrix 
model theory and topological string theory [26, 5]. In particular, one can obtain explicit formulae 
for the genus g free energies appearing in the 1/N expansion 



F{\gs) = Y^gl^-^F,{\), 

9=0 

where the 't Hooft coupling A is defined in (1.3), and 



(2.2) 



9s 



27ri 



(2.3) 



The genus g free energies Fg{X) obtained in this way are exact interpolating functions, and they 
can be studied in various regimes of the 't Hooft coupling. When A ^ they reproduce the 
perturbation theory of the matrix model around the Gaussian point. They can be also studied 
in the strong coupling regime A ^ oo, where one can make contact with the AdS dual. In this 
regime it is more convenient to use the shifted variable 



A = A 



1 

24' 



(2.4) 



As explained in [5], this shift is expected from type IIA and M-theory arguments [27, 28]. It 
turns out that, when expanded at strong coupling, the genus g free energies have the structure 



F,(A) = FP(A) + F7(A). 



(2.5) 



The first term represents the perturbative contribution in a', while the second term is non- 
perturbative in a', 

F7(A) ~ O fe-2^^^ (2.6) 
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and it was interpreted in [5] as the contribution of worldsheet instantons in the type IIA dual. 
For Fq 1 , the perturbative part is of the form, 



Ff(A) = |^-^log 



(2.7) 



2V2A 



while for > 2 one has 



where 



(2.8) 
(2.9) 



3=0 

is a polynomial. 

Besides the non-perturbative effects in a', one can use the connection between the large-order 
behavior of perturbation theory and instantons to deduce the structure of non-perturbative effects 
in the string coupling constant. In [8] a detailed analysis showed that these effects would have 
the form 

(2.10) 



exp ^— fc7r\/2A 



at large A. These were interpreted as D2-branes wrapped around generalized Lagrangian cycles 
of the target geometry. We will refer to these non-perturbative effects as membrane instanton 
effects, since they can be interpreted as M2 instantons in M-theory [29] but they are invisible in 
ordinary string perturbation theory. 

2.2 The partition function as an Airy function 

It was shown in [9] that the genus expansion of the perturbative free energies can be resummed. 
In order to do that, one has to use the variable [8] 



A. 



24 3fe2 



rather than (2.4). If we define the perturbative partition function as 

■^ABJM = 

then 



9=0 



2g-2 
s 



-^ABJM ^ 



A. 



where Ai is the Airy function. This can be also written in terms of N as 



where 



C-V3 



c 



k 



1 



24 3k 



TT^k' 



(2.11) 

(2.12) 
(2.13) 

(2.14) 

(2.15) 
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As noticed in [8], the expansion resummed in (2.14) makes perfect sense for finite k. Therefore, 
even if (2.14) was obtained from a calculation in the 't Hooft expansion, it should be part of the 
M-theory answer. Indeed, one of our goals in this paper is to verify this by computing 2^abjm 
directly in the M-theory expansion. 

The Airy function appearing in (2.14) gives an exact resummation of the long-distance ex- 
pansion in M-theory. To see this, one has to use the dictionary relating gauge theory quantities 
to gravity quantities. In particular, one has to take into account the anomalous shifts relat- 
ing the rank of the gauge group N to the Maxwell charge Q, which in turn determines the 
compactification radius L [27, 28]. The relation is 

« = "-2i(*-^)- 

The charge Q determines the compactification radius in M-theory according to 

L 



= 32Tr^Qk, (2.17) 

p/ 

where ip is the Planck length. The shift (2.11) was interpreted in [8] as a renormalization of the 
expansion parameter £p/L, since it means that the natural variable is 

- ^^/^ (2.18) 



L 



1 - 127r2 (ip/L) 



6 V6' 



and then the argument of the Airy function (2.14) is given by 

{256k.y^' (^i-^ . (2.19) 

The 1/N expansion of the ABJM matrix model was derived in [5] by using the holomorphic 
anomaly equations [30] of topological string theory. The result (2.14) was obtained in [9] by 
looking at the recursive structure of these equations. There is however a much simpler method 
to obtain (2.14) which exploits the wavefunction behavior of the topological string partition 
function^. Our derivation of (2.14) in this paper does not depend at all on ideas from topological 
string theory, but since it is formally very similar, we will now present this simpler argument. 
We will rely on results and notations of [5]. Readers who are not familiar with topological string 
theory can skip the rest of this section and proceed to the next one. 

As shown in [31], it follows from the holomorphic anomaly equations that the topological 
string partition function is a wavefunction on moduli space. In particular, its transformation 
from one symplcctic frame to the other is given by a Fourier transform. This property was 
spelled out in detail and exploited in [32]. The main result is summarized as follows. Let 



a (3 
7 6 



gSL(2,Z) (2.20) 



be a symplectic transformation relating two different frames (we assume for simplicity that there 
is a single modulus in the problem). This means that the periods {daFo, a) transform as 



oF I la 



(2.21) 



^We would like to thank C. Vafa for reminding us this. 
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Then, the full topological string partition function 



Z{a) = exp 



9=0 



transforms as 



where 



S{a, 



1 



-a7 



(2.22) 



(2.23) 



(2.24) 



In the context of ABJM theory, as explained in detail in [26, 5], the relevant quantities 
correspond to topological string computations in the so-called orbifold frame, where the natural 
periods are A (the 't Hooft coupling of the gauge theory) and the derivative dxFo. On the other 
hand, the most familiar frame in topological string theory is the large radius or Gromov-Witten 
frame, where the natural periods are T (the Kahler modulus) and the derivative OtF^^ ■ The 
genus g free energies in the large radius frame are given by the standard formulae. 



-^0 



fe>0 



kT 



fc>0 

k>0 

where Ng^k are Gromov-Witten invariants in the local x geometry (there is no constant 
term contribution at higher genus). The fact that the total free energy is at most cubic in T, up 
to exponentially small corrections, is a well known fact in topological string theory. 

In [5], the periods in the orbifold frame were written in terms of periods in the large radius 
frame in order to perform analytic continuations to strong coupling. By general principles, this 
relation must be a symplectic transformation like (2.20). In fact, it is easy to see that the results 
of [5] relating the periods can be written as the following symplectic transformation: 



-kT 



(2.25) 




where 



and 



A 



A= A, T=-T, 

c 2c 



27ri, 



Fo = Fo- ttHX, 



(2.26) 



(2.27) 



(2.28) 



Then, according to (2.23), (2.24), the total partition functions are related by the following for- 
mula: 

exp [F{X) - TT^iX/g^,] oc J dTexp [-f V^' + TX/g^, + ^^^(r)] . (2.29) 
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Notice that, up to nonperturbative terms in T, this is the integral of the exponential of a cubic 
polynomial, therefore we will indeed get an Airy function. Let us introduce the new variable n 
through 

T = ^ - TTi. (2.30) 

Then, one finds the expression 

/( 2u? k 1 / ' 



(2.31) 



oc 



Ai [(7-V3(Ar - B)] (l + C(e-2^^)) , 
where we used the following integral representation of the Airy function, 

Ai(^) = ^ exp (^j - z?j , (2.32) 

and C is a contour in the complex plane from e~^'^/^oo to e"^/^oo. In (2.31), C is given in (2.15) 
and 

The result of (2.31) is of course the expression obtained in (2.14). Notice that the first term in 
the shift B comes from Fq^^ in (2.28), while the second term is due to the first, perturbative 
term in . The exponentially small corrections in TV in (2.31), which are due to the world- 
sheet instantons at large radius of the topological string, become, after Fourier transform, the 
worldshcct instantons (2.6) of the type IIA superstring. 

This derivation is nice, but it seems difficult to generalize it in its current form to other 
Chern-Simons-matter theories, and prove in this way the conjecture (1.2) for other cases. In 
this paper we will find a completely different approach to the derivation of the Airy function 
which turns out to formally equivalent to the one based on topological string theory. However, 
this approach can be extended to many M > 3 CSM theories and makes it possible to verify the 
conjecture 1.1 for many of them. 



3. Chern-Simons-matter theories as Fermi gases 



3.1 ABJM theory as a Fermi gas 

Our Fermi gas approach is based on the following observation. The interaction term between the 
eigenvalues in (2.1) can be written in a different way by using the Cauchy identity: 



n,<,.[2sinh(^)] [2sinh(^^) 



n,,2cosh(^) 



= det. 



2 cosh (^) 



2 cosh 



(3.1) 



In this equation, Sn is the permutation group of N elements, and e(cr) is the signature of the 
permutation a. This identity has been used in other matrix models in [21, 22, 23] in order to 
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study them in the grand canonical ensemble, as we will do here. In the context of ABJM theory, 
it was used in [33] in order to prove the equivalence of (2.1) and the matrix integral for = 8 
super Yang-Mills theory in three dimensions, when k = 1. The manipulations in [33] can be 
easily generalized to arbitrary k, and one obtains the following expression for the ABJM matrix 
model, 

^^"^^ ^ ^ Sy^^'^' ^ (2-^)" n. 2 cosh (^) 2 cosh • ^'-'^ 



2k 



We will derive this expression below with a different technique, which can be used for more 
general Chern-Simons-matter theories. The main property of (3.2) is that it makes contact with 
the standard formalism to study partition functions of ideal Fermi gases. Indeed, let us introduce 
the function ^ ^ ^ ^ 

^^^^'^^^ = 2^ (2 eosh f ) (2 cosh ff"^ cosh (^ " ^^"^^ 
If we interpret it as a one-particle density matrix in the position representation 

p{xi,X2) = {xi\p\x2), (3.4) 

the matrix integral (2.1) can be written as the partition function of an ideal Fermi gas with A'^ 
particles 

^(^) = M E (-1)'^'^^/ d^xnp(^.^a«)- (3.5) 

It is well-known that the sum over permutations appearing in the canonical free energy of an 
ideal quantum gas can be written as a sum over conjugacy classes of the permutation group (see 
for example [34]). A conjugacy class is specified by a set of integers {m^} satisfying 

^imt = N. (3.6) 

e. 

Let us define 

= J dxi • • • dxi p{xi,X2)pix2, xs)--- p{xi^i,X(,)p{xe, xi). (3.7) 
Then, the partition function is given by, 

^W-E'n^^ (3.8) 

{ma t 

where the ' means that we only sum over the integers satisfying the constraint (3.6). 

Due to the constrained sum, the canonical partition function is not easy to handle for large 
N . As usual, the remedy is to consider the grand partition function 

oo 

S = 1 + ^ Z(iV)z^, (3.9) 



where 



(3.10) 
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Figure 1: A one-dimensional periodic lattice with l sites. The transfer matrix p can be regarded as the 
quantum propagator for a single particle in Euclidean, discretized time with Hamiltonian H . 

plays the role of the fugacity and /i is the chemical potential. The grand-canonical potential is 

J(^)=logH. (3.11) 

Notice that this potential (like the free energy) has the opposite sign to the usual conventions 
in Statistical Mechanics. A standard argument (presented for example in [34]) tells us that the 
sum over conjugacy classes in (3.8) can be written as 



J(/.) = -5^Z,^^. (3.12) 



{-zY 

The canonical partition function is recovered from the grand-canonical potential as 



At large N, this integral can be computed by applying the saddle-point method to 

Z{N) = ^J^f^ exp [J{fi) - ^lN] . (3.14) 

The saddle point occurs at 

^ = |^ = -E^^(--)^ (3-15) 

^ i>i 

and defines a function fi^:{N). The free energy is given, at leading order as — )• oo, by 

F{N) = J{fi^) - ^x^N. (3.16) 

However, it is possible to compute the 1 /N corrections to this relation by simply computing the 
corrections to the full integral in (3.14). This is what we will eventually do. Notice the similarity 
between the traditional inverse transform (3.14) and the Fourier transform (2.31) in topological 
string theory. 

We have then shown that the original ABJM matrix integral can be computed as the canon- 
ical partition function of a system of N non-interacting fermions, where the one-particle density 
matrix is given by (3.3). We just have to solve the corresponding one-body problem in order to 
compute the relevant thermodynamic quantities of the system. Equivalently, one should compute 
the quantity Zi introduced in (3.7). This quantity can be regarded as the partition function of 
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a classical lattice gas with I particles in a periodic lattice with nearest-neighbour interactions, 
as shown in Fig. 1. The density matrix p plays the role of the classical transfer matrix of the 
system (see for example chapter 12 of [35]). It defines a symmetric kernel 

(x|p|(^)= y" da;'p(a:,a;')'/'(a^'), (3-17) 

so that 

Z^ = Tr/. (3.18) 

It is easy to see that this kernel is a non-negative Hilbert-Schmidt operator, therefore it has a 
discrete, positive spectrum 

Pl^^n) = An|0n), W = 0, 1, • • • , (3.19) 

where are orthonormal eigenfunctions and we assume that 

Ao > Ai > A2 > • • • . (3.20) 
We can then write the density matrix as 

p = ^\n\<^n){<^n\- (3.21) 
n>0 

In terms of these eigenvalues we have, 

Z, = Y,\l (3.22) 

n>0 

When i is large, this sum is dominated by the largest eigenvalue Aq, 

Ag, £ » 1. (3.23) 

It also follows from this representation that the grand-canonical partition function is given by a 
Predholm determinant, 

S = det (1 + zp) = JJ (1 + z\n) . (3.24) 

n>0 

Instead of using the formulation of the lattice problem in terms of the density matrix oper- 
ator, we can introduce a quantum Hamiltonian in the standard way, 

p = e-". (3.25) 

This leads to the well-known equivalence between the partition function of a classical lattice gas 
(3.18) and the propagator of a quantum particle in i units of discretized time (see for example 
[35, 36]). We can then write 

Z^ = Tre-^^. (3.26) 

To find the Hamiltonian corresponding to the ABJM matrix model, we first write the density 
matrix (3.3) as 
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In this equation, q^p are canonically conjugate operators, 

[q,p] = ih, 

and 

h = 27rA;. 



(3.28) 



(3.29) 



This is a key aspect of this formahsm: h is the inverse coupling constant of the gauge theory/string 
theory, therefore scmiclassical or WKB expansions in the Fermi gas correspond to strong couphng 
expansions in gauge theory/string theory. The potential U{q) in (3.27) is given by 



C/(g) = log(2cosh|), 



and the kinetic term T(p) is given by the same function. 



T{p) = log ( 



2 cosh ^) . 
2/ 



(3.30) 



(3.31) 



The peculiar kinetic term (3.31) can be regarded as a non-trivial dispersion relation interpolating 
between the quadratic behavior of a non-relativistic particle at small p, 



T(;>)-log(2) + |-, p^O, 
and the linear behavior of an ultra-relativistic particle at large p, 

T{p) ~ ^, \p\ 



GO. 



(3.32) 



(3.33) 



Notice that, as it is standard for Hamiltonians defined by transfer matrices at finite lattice 
spacing [35, 36], the quantum operator H defined by (3.25) and (3.27) differs from 



T{p) + U{q) 



(3.34) 



in h corrections. There is a very elegant method to obtain these corrections based on the phase- 
space or Wigncr approach to quantization. This method will be also extremely useful in setting 
the semiclassical or WKB expansion of our thermodynamic problem. We first recall that the 
Wigner transform of an operator A is given by (see [37] for a detailed exposition of phase-space 
quantization) 



^w(g,p) 



dq Iq 



A 



q+-)e 



(3.35) 



The Wigner transform of a product is given by the ^-product of their Wigner transforms. 



AB = * Bw 

w 



where the star operator is given as usual by 

★ = exp 



ih 

2^ 



and is invariant under linear canonical transformations. Another useful property is that 

dpdq 



TvA 



27rh 



Aw{q,p). 



(3.36) 



(3.37) 



(3.38) 
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In order to calculate the h corrections to the Hamiltonian, we consider the Wigner transform of 
the density matrix (3.27). By using (3.36) we find, 

pw = e-^^(«)*e-^(f)*e-^^(«). (3.39) 

Let us note that the partition function depends only on the eigenvalues A„ of p (or, equivalently, 
on the traces = Tr p^). Therefore there is the following freedom in the choice of p: 

p VpV-^ (3.40) 

which translates into 

Pw{q,p) yw{q,p) -k pw{q,p) -k {V~^)w{q,p) (3.41) 

after the Wigner transform. Equation (3.39) defines the Wigner transform of our Hamiltonian 
through 

/>w = e-^w, (3.42) 

where the ★-exponential is defined by 

e^p^{A) = l + A+^A-kA + ---. (3.43) 

The quantum Hamiltonian iJw can be computed by using the Baker-Campbell-Hausdorff for- 
mula, as applied to the ^-product. One finds, 

HMq,p) =T + U+^[T,[T, UH + ^ [U, [T, UH + ... 

= T{p) + U{q) - ^ {T'{p)f U"{q) + |^ {U'{q)fT"{p) + 0{h\ 

where we have used the fact that, at leading order in h, the Moyal bracket is the Poisson bracket 

[A,B]^ = Ai<B-Bi.A = ih{A, B} + 0(h^). (3.45) 

Further corrections to (3.44) can be computed to any desired order, see (A.l) for the result at 
order 0{h^). 

3.2 More general Chern Simons matter theories 

The identification of the matrix model of ABJM theory as the partition function of a Fermi gas 
can be also made for more general A/" > 3 Chern-Simons-matter theories. We will set up the 
formalism for the necklace quivers with r nodes considered in [38, 39], and with fundamental 
matter in each node (see Fig. 2). These theories are given by a 

U{N)k, X U{N)k^ X • • • U{N)k^ (3.46) 

Chern-Simons quiver. Each node will be labelled with the letter a = I,-- - ,r. There are 
bifundamental chiral superfields Aaa+i, ^oa-i connecting adjacent nodes, and in addition we 
will suppose that there are Nf^ matter superfields {Qa, Qa) in each node, in the fundamental 
representation. We will write 

ka = riak, (3.47) 
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Figure 2: A quiver with, r nodes forming a ned?;lace. 



and we will assume that 



a=l 



(3.48) 



According to the general localization computation in [2], the matrix model computing the 
partition function of a necklace quiver is given by 



Z{N) 



m 



nyJU 



dA„.,; exp 



47r ^a,i 



2-K 



n 



n 



i<j 



2 sinh 



(2 cosh ^i)"^'" ail 2 cosh (^^^^^^) ■ 



(3.49) 



This matrix model is very similar to the Ar-i models considered in for example [22, 40], and one 
can use a very similar strategy in order to rewrite them as Fermi gases. First of all, we define a 
kernel corresponding to a pair of connected nodes (a, b) by, 



Kab{x',x) 



1 ^^P{il^} 

27rfe2cosh(3^) 



2 cosh 



2k 



-N 



fb 



(3.50) 



where wc set x = X/k. The grand canonical partition function corresponding to the above matrix 
model is defined as in (3.9). Then, if we use the Cauchy identity (3.1), a simple generalization of 
the above arguments makes it possible to write it again as a Fredholm determinant (3.24), where 
now [22] 

p = KrlKi2 ■ ■ ■ Kr-l,r (3.51) 

is the product of the kernels (3.50) around the quiver. Therefore, we can apply exactly the same 
techniques that we used before in ABJM theory. In a sense, we are "integrating out" r — 1 nodes 
of the quiver in order to define an effective theory in the r-th node, but with a complicated 
Hamiltonian which takes into account the other nodes. 
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This idea can be made very concrete by looking at the Wigner transform of the operator p 
in (3.51). We first compute the Wigner transform of the kernel (3.50), 

where the h in the ★ product is given again by (3.29). Let us note that 

ad* 



e 2a ★/(p)*e 2^ = / 2ft -kpi^e 2fi J = / L ^pj=f{p — nq), (3.53) 
where we used that 

[q^,p]^ = mq. (3.54) 
We obtain then, for the Wigner transform of the density operator (3.51) 

1 1 1 



pw{q,p) 



2 cosh I [2 cosh A] 2 cosh M 
1 1 1 

[2 cosh A] ^^'2 * 2cosh ^-(%+"^)^ " * [2 cosh 4] ^^3* (^•5^) 

1 1 

2 cosh Pr(ni+L -^+nr-r)g * ^^^^^ jV/. 

where we used (3.48). For necklace theories without fundamental matter this is simply 

1 1 1 1 

^^^'^^ P> - 2^ ^ 2 cosh ™ 2 cosh '"'"'2 cosh P^-+n.-^)i • ^^-^^^ 

In particular, for the ABJM necklace {—k, k) with fundamental matter Nf^ = Nf^ = Nf first 
considered in [41, 42, 43], we have 

« ^' = 2^ * [2 , J ^f,* i^SPH * [2 cosh if 

If we perform a canonical transformation 

P ^ —1, q—^p + q (3.58) 

and we conjugate by [2 cosh |]^^^ to obtain a symmetric kernel, we get the equivalent represen- 
tation, 

"^^'"'^ ^ [2cosli]^/^ ' [2cosh'l^]^^ * 20^* [2cosh'^]^^ " [2cosli]^/^ ^'"''^ 

which, for Nf = 0, agrees with the result (3.39). In this way, we have reduced the general 
necklace quiver theory to an ideal Fermi gas whose one-particle quantum Hamiltonian is defined 
by the above density matrices through (3.42). 

Notice that, in general, the density operators p are not Hermitian, and correspondingly 
is generally not real. This reflects the fact that the free energy on the three-sphere of these CSM 
theories is in general complex. 



- 16 - 



4. Thermodynamic limit 



It is well-known that the thermodynamic limit of an ideal quantum gas can be evaluated by 
treating the one-particle problem in the semiclassical or WKB approximation. Moreover, the 
corrections to the thermodynamic limit can be obtained by studying the quantum corrections to 
the semiclassical limit. In this section we will present general results about the thermodynamic 
limit and we will illustrate them in ABJM theory. More general theories will be considered in 
section 6. 

4.1 The thermodynamic limit of ideal Fermi gases 

In the following we will need several standard results in the analysis of ideal quantum gases. The 
distribution operator at zero temperature is given by, 

h{E) = e{E-H) (4.1) 

where 9{x) is the Hcaviside step function. The trace of this operator gives the function n{E), 
counting the number of eigenstates whose energy is less than E: 

n{E) = Trh{E) = J2^{E-En)- (4.2) 

n 

Notice that 

£;„ = -logA„, (4.3) 
where A„ are the eigenvalues (3.19) of the density matrix. The density of eigenstates is defined 

by 

p{E) = '^ = Y.5iE-Er,). (4.4) 

n 

The one-particle canonical partition function is then given by the standard formula, 

/•oo 

Zi= / dEp{E)e-^'', (4.5) 
Jo 

while the grand- canonical potential of the N particle system is given by 

/•oo 

J(m)= / dE p{E) log [1 + ze-^) . (4.6) 

Let us now consider the thermodynamic limit of the system, when N ^ oo. In this regime, 
the behavior of the system is semiclassical and the spectrum of the one-particle Hamiltonian is 
encoded in the functions n{E), p{E). The thermodynamic limit is governed by the behavior of 
these functions as £^ S> 1. We notice that, if 

n{E) ^ CE', E » 1, (4.7) 

then the grand-canonical potential is given by 

oo 

J{p) ^sC J log (1 + ze-^) E'-^dE = -CT{s + 1) Li^+i(-e''), (4.8) 
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where Li^ is the usual polylogarithm function. The number of particles is related to the chemical 
potential by 

N{n) ^ Cr{s + 1) Li,(-e'*), (4.9) 
and large N corresponds to large /j,. In this regime, we have 

JifJ-) ~ -^t^'''^ ^(/^) ~ Cf^' ■ (4-10) 

S ~r -L 

The second equation defines n as function of A'^, and we deduce from (3.16) that the canonical 
free energy is given, as TV — >■ oo, by 

F{N) « -^^<^"^^'^^ • (4.11) 

These formulae should be familiar from the elementary theory of ideal quantum gases. For 
example, the textbook ideal Fermi gas in three dimensions has s = 3/2. 

To determine the value of s for a given system we notice that, in the semiclassical limit, the 
trace is replaced by an integral over phase space 

Tr^f^ (4.12) 

which gives the standard semiclassical formula 

i.e. the number of eigenstates is just given by the volume of phase space. The surface 

H{q,p)=E (4.14) 

is just the Fermi surface of the system. For a one-dimensional ideal gas whose one-particle 
Hamiltonian is of the form 

H r-. AIpI'' + B\qf (4.15) 

we have 

.= 1 + 1. (4.16) 

This will be useful later on. 

4.2 A simple derivation of the N^/^ behavior in ABJM theory 

We can now study the thermodynamic limit of the partition function of ABJM theory. In this 
case, the Hamiltonian appearing in the semiclassical formula (4.13) is just given by the classical 
counterpart of (3.34), 

H^l{q,p) = Tip) + U{q) = log (2 cosh |) + log (2 cosh |) . (4.17) 

Here we have neglected the h corrections appearing in i?w It is easy to show that the minimum 
energy is 

^0 = 2 log 2, (4.18) 
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Figure 3: The Fermi surface (4.20) for ABJM theory in the q-p plane, for i? = 4 (left) and E — 100 
(right). When the energy is large, the Fermi surface approaches the polygon (4.22). 

which corresponds to the maximal eigenvalue of the density matrix 

Ao = \. (4.19) 

This is the semiclassical value given by the leading WKB approximation, and it will be corrected 
quantum-mechanically. In the large E regime, the discrete spectrum "condenses" along a cut in 
the complex plane, and Ao signals the endpoint of the cut. 

In order to proceed with the analysis of the thermodynamic limit, we should determine the 
Fermi surface 

Hd{q,p)=E (4.20) 

controling the density of eigenvalues. We show the shape of this surface in Fig. 3 for E = 4 (left) 
and E = 100 (right). It is clear that in the thermodynamic limit, when E is large, the surface 
can be approximated by considering the values of U{q), T{p) for q,p large. In this regime we 
have 

U{q)^^-, kl^oo, T(p)«M, N^oo, (4.21) 
so that (4.20) is approximately given by 

\q\ + \p\ = 2E, (4.22) 
as it is manifest in the graphic on the right in Fig. 3. From (4.15) and (4.16) we deduce that 

s = 2. (4.23) 



Since 

the number of states is given by 
By comparing with (4.7), we find 



Vol(^) w (4.24) 
niE) « -^E\ (4.25) 
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The equation (4.11) gives immediately 



F(iv)«.-^^^iV3/2. (4.27) 

This is exactly the result first found in [5] using the 't Hooft expansion of the matrix model. 

The derivation presented here is however completely elementary, and relies on basic notions of 
quantum Statistical Mechanics: the 3/2 scaling of the number of degrees of freedom is nothing but 
the scaling of the free energy of an ultrarelativistic gas of one-dimensional fermions in a linearly 
confining potential. No matrix model techniques are needed. In this sense, our derivation is even 
simpler than the one presented in [10], which required some detailed analysis of the eigenvalue 
interaction in the matrix integral. 

We would like to emphasize that the above result (4.27) provides the right large N behavior 
of the system at finite k. This is because the true expansion parameter in the semiclassical 
expansion is h/E, which is small for large E even at finite h. This can be proved rigorously for 
some spectral problems defined by kernels of the form (3.3) [44], and we will verify it in section 
5 by a detailed analysis of the WKB expansion. 



4.3 Large N corrections 

One advantage of the statistical-mechanical framework presented here is that it makes it possible 
to compute corrections to the thermodynamic limit in a systematic way. To start the study of 
these corrections, we now look at the thermodynamics of the Fermi gas of ABJM theory in the 
semiclassical approximation, but taking into account the exact value of the volume of phase space 
(i.e. we go beyond the polygonal approximation in (4.22)). As expected, this gives sub-leading 
and exponentially suppressed corrections at large A^. 

The computation of the exact volume is equivalent to computing all the exactly in the 
semiclassical approximation, and resumming the resulting series (3.12). Using that 

^2cosh|j ^^^^ 



we find 



where 



Therefore, 



where 



Zi^-Z'^''', (4.29) 

(0) _ f dpdg _ 1 r\e/2) 



J{li) « ^ Jo(m) (4.31) 



Jo(/.) = -£ 



^^^7f2(£)" 
1 „/lll 2> ^2 „/_33 



4^^^M2'2'2'-^'2' 16 J 87r2^^^ r ' ^' ^' ^' 2 ' 2 ' 16 



(4.32) 
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This function has a branch cut in the 2;-plane at (—00, —4]. This is expected: indeed, from (3.24) 
there should be a cut starting at 

z = -Ag ^ = -4, (4.33) 

indicating the condensation of eigenvalues for the one-particle density matrix. The function 
(4.32) has the following asymptotics for large /x, 

The leading, cubic term in is the responsible for the behavior (4.27). The sublcading term in n 
gives a correction of order N^^"^ to the leading behavior (4.27). The last, non-perturbative term 
involves an infinite power series of exponentially small corrections in /x. They have the structure, 

00 

e=i 

Explicitly, one finds for the very first orders 

J^P(^) = ^ (6 - TT^ + 6/x - 6(1^) e-^f" + 7^ (25 - Gtt^ - - 36/i^) e'"^" 

3tt ztt (4.36) 

+ O ifi^e-^^") . 

The non-perturbative part of the grand potential leads to exponentially small corrections in N in 
the canonical free energy. In fact, using (4.10) we find that, once evaluated at the saddle-point, 

exp (-2/x) « exp [-sPlT^k^l'^N^'^^ . (4.37) 

This is precisely the action for membrane instantons (2.10) found in [8] as large N instantons of 
the matrix model in the 't Hooft expansion. We conclude that the exponentially small corrections 
in which in this approach appear already in the semi-classical approximation, correspond in 
fact to non-perturbative corrections in the genus expansion, and should be identified as membrane 
instanton contributions. 

As mentioned before, the calculation of these exponentially small corrections to the grand- 
canonical potential is equivalent to the exact calculation of the volume (4.13) of classical phase 
space. To see this, we notice that we can write this volume as a period of the one-form pdq along 
the curve (4.14) 

Vol(E) = j> pdq. (4.38) 

H^x{q,p)=E 

This period vanishes at the point E = Eq. It turns out that its exact value is given by a Meijer 
function, 

E / JIE 1 1 1 \ 

This leads to the following large E expansion of the number of states, 

2 1 

n{E) = — --+ 0{Ee-'^) + 0{h), (4.40) 

where the first term agrees of course with the semiclassical calculation at large E done before. 
One can then check that the expression (4.6) for the grand-canonical potential reproduces (4.32), 
once the density obtained from (4.39) is used. 
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4.4 Relation to previous results 

The semiclassical limit of the one-particle Hamiltonian turns out to be closely related to the 
planar limit of ABJM theory studied in [26, 5, 8]. 

First of all, the semiclassical quantization of the one-dimensional problem leads to the Fermi 
surface 

T{p) + U{q) = E (4.41) 
which is in fact a curve in phase space. Let us now make the following change of variables, 

a; = | + |, y = p + Tri, (4.42) 
which, up to an overall constant, preserves the form pdq. In terms of the exponentiated variables 

X = ei/2+P/2^ Y = -ef . (4.43) 

the Fermi surface (4.41) reads 

y + — - + IkX - 1 = 0, (4.44) 

where 

IK = e^. (4.45) 

The curve (4.44) is nothing but the spectral curve of the ABJM matrix model written down 
in for example [5]. The minimal energy Eq given in (4.18) corresponds to the conifold point 
K = 4i studied in detail in [5]. The volume of phase space, which as we remarked after (4.38) 
is a vanishing period at E = Eq, is actually proportional to the conifold period studied in [8]. 
Finally, the large energy limit, in which the Fermi surface becomes a polygon, is nothing but the 
tropical limit of the spectral curve, studied in [11]. 

5. Quantum corrections 

In the previous section we have recalled the semiclassical limit of ideal Fermi gases, and we have 
studied in detail the case of ABJM theory. We now study the corrections to the semiclassical 
limit in a systematic and general way. These corrections lead to a power series in h? oc k'^ for 
the grand- canonical partition function. As we will see, only the first correction contributes 
to the asymptotic series in 1/N of the canonical free energy, up to an additive function of k but 
independent of N. This means that we can compute the full series of 1/N corrections to the 
original matrix model partition function, up to an overall, A^-indcpendent constant. However, 
the exponentially small terms in appearing in J(/x) receive corrections to all orders in h?. 

5.1 Quantum-corrected Hamiltonian and Wigner— Kirkwood expansion 

There are two sources of h corrections in the one-body problem appearing in our Fermi gas for- 
mulation. The first one appears already in the Hamiltonian H: when we compute H starting 
from (3.27), the non-commutativity of the operators in (3.27) leads to 0{h) corrections to (3.34). 
This first source of corrections is nicely encoded in the Wigner transform (3.44). Another source 
of corrections is due to the standard semiclassical expansion of the density of eigenvalues. We 
now present a formalism to treat in a systematic way both types of corrections. This formalism 
is a generalization of the standard Wigner-Kirkwood h expansion [45, 46] in quantum statistical 
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mechanics, and it incorporates general, ?i-depcndent Hamiltonians. As in section 3.1, the for- 
mahsm is most conveniently formulated in the phase-space approach to quantization, and it has 
been developed in the context of many-body physics. The most elegant presentation is due to 
Voros [47, 48] (see also [49]). 

Let H be the Hamiltonian of a one-particle, one-dimensional quantum system, and let Hyj be 
its Wigner transform. We would like to compute systematically the h expansion of the canonical 
partition function and of the density of states. Following [48] we notice that it is possible to 
expand any function f{H) of H around HT^{q,p), which is a c-number. This gives, 



r>0 



(5.1) 



The semiclassical expansion of this object is obtained simply by evaluating its Wigner transform, 
and we obtain 

/(^)w = E V^'^ (^w)^. 



r>0 



where 



w 



and the Wigner transform is evaluated at the same point q,p. Of course, one has 

go = 1, Qi = 0, 



(5.2) 
(5.3) 
(5.4) 



and the quantities Qr for r > 2 can be computed again by using (3.36). They have an h expansion 
of the form 

gr= ^^"^^^"^ ^>2- (5.5) 

">[^] 

This means, in particular, that to any order in fi^, only a finite number of QrS are involved. One 
finds, for the very first orders [48, 49], 



92 = - 

93 = - 



4 



w 



gq2 Qp2 \y dqdp ) 

2 pa , 



/dHwy (THw^ f dHw \ ^ d^Hw ^ dH^ dHw d^H^ 
\ dq J dp^ \ dp J dq^ dq dp dqdp 



(5.6) 



+ 0(^4). 



One can then apply this method to compute the semiclassical expansion of any function of the 
Hamiltonian operator. For example, when applied to (4.2), one finds, 



r=2 ^' 



(5.7) 



therefore 



n{E) = f 
Jh 



dqdp 



+i:^/||e^*<'-*(^-«w). (5.8) 



'Hw(q,p)<E 27rn 

When applied to the canonical density matrix at inverse temperature /3, one finds. 



(5.9) 
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The standard Wigner-Kirkwood h expansion of the canonical partition function [45, 46] is just 
a particular case of (5.9) when the Hamiltonian is 



(5.10) 



Let us now apply this formalism to our case. First of all, the quantum-corrected Hamiltonian 
is given by a power series in h of the form. 



2n u{n) 

W ■ 



n>0 



At leading order we find of course the classical Hamiltonian (4.17), 

h'^ =T{p) + U{q\ 



(5.11) 



(5.12) 



the 0{h?) term is written down in (3.44), and the O^hf^) term can be found in (A.l). The 
one-particle canonical partition function can then be computed as a power series in h, 



-. oo 



n=0 



where 



.(0) 



dqdp 



(5.13) 



(5.14) 



is the classical limit. The expansion is obtained by grouping fi^ corrections in the expression 



The power series in h for Z(, leads to the following power series in A; for J(/x), 

^ oo 
n=0 

where 



(5.15) 



(5.16) 



(5.17) 



As an illustration of the above, general considerations, we will now calculate the first, 
correction to the semiclassical result of ABJM theory obtained in section 4.3. Using the formulae 
above, we find 



-IK 



(1) 
w 



2 6 • 



(1) 



-in,: 



+ 



dgdp 



24 



l(t/'(g))2r"(p)-l(r'(p))2C/"(g) 



{(Jj\q)fT"{p) + Tj'\q){:r\v)f\ - -U"{q)T"{p) 



(5.18) 
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To evaluate these coefficients, we need the integral appearing in (4.28), as well as 

tanh^(e/2) _ r^{i/2 + 1) 

7-0O ^2cosh(e/2))^ r(^) r{e + 2) ' ^'-^'^ 



We then find, 

T2(£/2 + l)r2(£/2) T\£/2 + l)] e T\t/2 + l) 



4r(^ + 2)r(^) r2(^ + 2) 



167r r2(^ + 2) 



(5.20) 



Prom (5.20) one can compute Ji{z) in closed form. Let us introduce the function 

/ 3 3 z^X 2^ / 3 5 

/(.) = 3i^2 ( 1, 1, 1; -, -; - j - ^ 3F2 h , 1, 2; -, -; - 



(5.21) 



+ 1 (-2.^(1)-. + .^) 
where E{k) is the complete elliptic integral of the second kind with modulus k. Then, one finds 

^i(/^) = ^|/(^)-(^^)'/(^)|- (5.22) 
The asymptotic expansion of the above function at large /i is given by 

f{z) - (^^) ' /(^) = M - 2 + O (M'e-2^) . (5.23) 

Therefore, we find, at next-to-leading order in A;, the following expression for the grand canonical 
potential of ABJM theory. 

Notice that the non-perturbative corrections in |Lt to (5.23) involve only even powers of z. This 
is consistent with their interpretation as membrane instantons. 

5.2 General structure of quantum corrections 

As we mentioned above, we can compute the quantum corrections to J(//) either by working out 
the corrections to the Zi integrals, or by working out the corrections to the function n{E). In 
order to understand the general structure of these corrections for a Fermi gas, the second point 
of view is more convenient. In this section we will analyze this general structure in detail, and 
we will make a precise connection between the structure of n{E) and the expected Airy function 
behavior. 

First of all, we have to understand more precisely the relationship between the structure of 
n{E) and the structure of J{ijl). Let us write the density function n{E) in the form, 

n{E) = CE'^ + nQ + n^^{E), (5.25) 

where the last term has the following asymptotics at infinity, 

nnp(-E) = 0(£;e--^), E ^ 00. (5.26) 
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We know from (4.40) that this is indeed the case at leading order in k for ABJM theory and in 

the next subsection we wih show that quantum corrections do not spoil this behavior. Notice 
that, since all eigenvalues of our Hamiltonian are positive, we must have 

n{0) = 0, (5.27) 

therefore 

n„p(0) = -no. (5.28) 

If we now plug (5.25) in (4.6) we find. 



J(^) = J dEn'{E) log(l + ze-^) 



oo oo oo 

= -2CLi3(-z) + M j dEn'^piE)- j dEn'^^{E)E + J dE n'^^{E) \og{l + e"" / z) . 



(5.29) 





The second integral gives 



/•oo 

/ dEn'iE) = -nnp(O) = no, (5.30) 
Jo 

where we used (5.28). The last term can be calculated as 

oo oo 

J dEn',^{E)logil + e^/z) = mlog{l + l/z)- J dE^"^^, (5.31) 



and both terms are non-perturbative in /x. Indeed, 



oo 

J^^^-J^^T^-0(,.-n. (5.32) 





Then, by using the standard asymptotics of the trilogarithm 



Li3(-^) = -^-yM + 0(e-'^), (5.33) 



we deduce the following asymptotic expansion of J(/x) for large /x: 

C 
3 

where 



Jil^) = ^^1^^ + Bfi + A + Jnpin) , (5.34) 



B = no + 



oo 

A = -Tr'H = -J dEEn'^^{E), 



(5.35) 



and 

Jnp(/x) = 0(Me-'^). (5.36) 
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Notice that A is a non-trivial function of k, but it doesn't depend on fj,. If we now plug this in 
(3.14), we find immediately 



Z{N) = C-^/^e^ Ai C-^/^{N - B) + Z^piN), 



(5.37) 



where the last term is non-perturbative in N. 

We then see that, if we are able to derive the structural results (5.25) and (5.26) for the 
density of states of a given theory, the conjecture 1.1 for the M-theory expansion is proved. In 
fact, so far we have not specified in which regime we are working in k. In practice, we have to 
work in an expansion in k around A; = 0. However, we expect that C will only get contributions 
at leading order in k (i.e. the strict scmiclassical limit), and that B will be only corrected 
at the ncxt-to-lcading order in k. Wc will now verify this in ABJM theory. In contrast, the 
//-independent term A gets corrected at all orders in k. 

5.3 Quantum corrections in ABJM theory 

We now study the general structure of quantum corrections in ABJM theory, by using the strategy 
explained above, i.e. by looking at the number of eigenvalues n{E). Our goal is to show that 
n{E) has the structure (5.25). This involves a somewhat detailed argument. Since not every 
reader might go through it, we want to emphasize that the physics behind this argument is very 
simple. The WKB expansion of the density of eigenvalues of a quantum system is in fact an 
expansion in 



Therefore, if the leading order term in n{E) is of the form CE^, the first quantum correction 
gives the constant term in (5.25), and further terms in the WKB expansion do not correct the 
polynomial part of n{E). They can only give exponentially small corrections in E. In the rest 
of this section, we will verify that this qualitative argument is actually correct in the case of the 
one-body problem appearing in ABJM theory. 

Our starting point in the study of quantum corrections in ABJM theory is (5.8). As we know, 
there are two sources of h corrections in this formula. One is the quantum-corrected Hamiltonian, 
and the other are the terms Qr appearing in the generalized Wigner-Kirkwood expansion. We 
will consider first the quantum corrections coming from Ht^, i.e. from the first term in (5.8). 
Since we have a symmetry q —q and p — >■ —p in the problem, we can restrict ourselves to the 
case q^ and p^ 0. We want to solve the equation 



in the limit E ^ oo. This defines a "quantum curve" or "quantum Fermi surface," including 
explicit h corrections. At leading order in E the curve is given by (4.22), and the corresponding 
domain (in the positive quadrant) has volume 




(5.38) 



Hw{q,p) = E 



(5.39) 



Volo(^) = 2E'^. 



(5.40) 
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One crucial ingredient in what follows is the fact that the function U {q) and its derivatives have 
the following asymptotics as q ^ 1: 



Q , (-1) 



[/(g) = log2cosh| = |+j; 



fc+i 



2 2 ^ k 



-e 



-kq 



U'{q) = \ tanh l = \ + ^(-1) V^'', (5.41) 



4 cosh"^ 2 
4COsn 2 

The same results hold for T{p). Notice that, if we take a number large enough of derivatives of 
these functions, they become exponentially suppressed at infinity. This will be eventually the 
source of the simplifications at large E. 




Figure 4: The regions I (left) and II (right) under the quantum curve H^{q^p) ^ E in the positive 
quadrant. The diagonal dashed line is the polygonal curve (4.22). 

Let us now consider the point {q*,p*) in the curve (5.39), where 

p* = E. (5.42) 

It is easy to see from the explicit form of i^w that 

q^=E + (e-^) (5.43) 

where the exponentially small corrections in E are themselves power series in h?. This point 
divides the curve (5.39) into two segments, and defines two regions for the fully corrected volume, 
as shown in Fig. 4. Region I is defined as the region under the quantum curve with p > p^,, while 
region II is defined by g > . We have 

Vol(^) = 4Voli(E) + 4Volii(S), (5.44) 

where 

Yoli{E)= / p{E,q)dq, Vol„(^) = / q{E,p)dp - p,q,, (5.45) 
Jo Jo 

and p{E,q) and q{E,p) are local solutions of H\^{q,p) = E. 
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Let us first consider the curve bounding region I. Along tliis curve, p{q, E) > E, tlierefore 
exponential terms in p in i^w are bounded by exponential terms in E. We can then write 

T{p) = I + O (e-^) , T'ip) = \ + (e-^) , (5.46) 

and 

T(")(^?) = C»(e-'^), n>2. (5.47) 

In the quantum corrections to the function we will have terms of the form {T^^^p))"^ , with 
k > 1. Due to (5.46) and (5.47), and neglecting exponentially small corrections of the form 
O (e~^), we should keep only the terms (T'(p))" with n > 1 (like the third term in the second 
line of (3.44)). But these terms always multiply terms of the form U^'^'^\q). We conclude that, 
on the curve bounding region I, 

Hw = l + U{q) - ^U"{q) + ^ ^ h''^CnU^''^\q) + O (e"^) . (5.48) 

n>l 

The third term in this expression comes from the third term in the second line of (3.44). The 
fourth term comes from higher quantum corrections (see the first term in the last line of (A. 2) 
for an example of such a term at order 0{h^)). We can now solve for p along this curve, 

p{E,q) = 2E-q + Ap{E,q), (5.49) 

where 

ApiE, q)=q- 2U{q) + ^^"(g) - ^ h^^CnU^^^\q) + ©(e"^). (5.50) 

n>l 

We calculate the volume of region I as follows, 

Voli = Vol? + AVoli. (5.51) 
The first term comes from the polygonal limit of the curve, 

Yo\f{E)= j {2E-q)dq = 2Eq,{E)-^^. (5.52) 


The second term comes from the corrections to the curve, and it is given by 



AVoli(E)= J Ap{E,q)dq 




-2 j (c/(g)-|)dg+^ j U"{q)M-Y,^^''cn j U('^\q)+0{Ee-^) 
"^^ 

= -2 J {U{q)-l) dq+^J U"{q)dq -^h'^Cn J U'^'^Hq) + O {Ee'^) 
"^-'^ 

(5.53) 
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In the last calculation we used that, up to non-perturbative terms in E, we can extend the 
integration region to infinity, and also that 



J U^'^'"'\q)dq = J7(2«-i)(oo) - U^'^'^-^\{S) = forn > 1. (5.54) 


A similar calculation can be done for region II. We obtain, from the polygonal approximation of 
the curve, 

\ol\i{E) = 2Ep,{E) - - p,{E)q,{E), (5.55) 



while the corrections give, 

AVol„(£;) = - ^ + O (Ee-^) . (5.56) 

Using that 

{E) + {E) = 2E + {e-^) (5.57) 

we finally get 

YoliE) = 8E^-^ + ^ + {Ee-^) . (5.58) 

We now consider the contribution from the quantum corrections to the density. In fact, these 
terms only give non-perturbative corrections in E. Using that 

SiE-HM,,p)) = S{p--p{E,q)) /^-^'^=5{q-qiE,p)) j ^J^^ (5.59) 

one can always decompose an integral over the phase space as a sum of one-dimcnsional integrals 
in regions I and II, as in (5.44). For region I one can use again the expression (5.48) and the 
properties (5.46), (5.47). The only nontrivial term which gives an correction comes from 
and gives. 



r dH^{q,p)d^H^{q,p) 



I 



dg 



24 5£;2 J ^ dp dq^ 




P=v{E,q) 



24 dE"^ 



[ dQ\Yl ^^''cnU^^''+'^\q) + 0(e-^) j = 0{Ee-'^). (5.60) 
\">o / 



For higher order corrections everything that contains a term with 

d''H^M{q,p) 



dp"" 

or with 



r > 1, (5.61) 
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is of order e ^ . Since the derivatives dp and dq always come in pairs, the only terms possibly 
contributing are of the form 

ill dp ) -ii g^. +^^^ ) ^^-^^^ 

i i 

where > 2, Vi After integrating and applying d^/dE'^ this gives a correction of order 0{Ee~^) 
by the same reason. 

We conclude that, to all orders in the h expansion. 

Therefore, by using (5.34) and (5.35) we find the expression 

Jabjm()u) = ^ + ^ (4 + ^) + ^('^^ + '^"P^'") ^^-^^^ 

where 

oo 

Jnp(M) = (a^,„/x' + fe^n/x + Q,n) fc^^-^^e'^^''. (5.66) 

e,n=i 

It is not manifest from the above results that this series involves only even powers of z~^, but 
we have verified it to be the case for the first three orders in k, and we believe it is a general 
feature. Finally, it follows from (5.64) and (5.37) that 



^ABJM(iV) = C-V3e^('=)Ai 



+ Znp(iV), (5.67) 



where C is given in (2.15) and Znp{N) are exponentially suppressed corrections at large A^. 
This concludes our derivation of the Airy behavior for ABJM theory. The function A[k) can in 
principle be determined, order by order in A;, by computing the Z^'^\ rcsumming the resulting 
series, and expanding at )U = oo, as we did in sections 4.3 and 5.1. One obtains, 

A sketch of the computation leading to the third term of this expansion can be found in the 
Appendix. In principle one can also compute A{k) by using the representation (5.35). What is 
the interpretation of A(A;)? One natural possibility is that A{k) encodes effects of order 

O (e-'^) - O (e-^/^^^ , (5.69) 

i.e. that A{k) gives the contribution from DO branes. Notice that the second and third coefficients 
in (5.68) are given by 

for g = 1,2. These are the coefficients of the power series expansion in k of 



2 r'' 



vr / ^2 log 
^ Jo s 



sm(e/2) 

e/2 



(5.71) 
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Figure 5: The Riemann surface of Hc\{p,q) = E for ABJM theory for large E. The four interior tubes 
form the limiting polygon of the Fermi surface. 

This gives indeed corrections of the form (5.69). It would be interesting to verify that the above 
function resums the small k expansion of A{k)^ . 

5.4 Quantum-mechanical instantons as vi^orldsheet instantons 

One obvious question that one can ask at this point is the following: where are the worldsheet 
instantons (2.6) that one finds perturbatively in gs in the 't Hooft expansion? We now give 
some preliminary evidence that worldsheet instantons correspond to the quantum-mechanical 
instantons of the Hamiltonian -ffw 

So far our focus has been in the perturbative corrections in h, but one should expect gener- 
ically non-perturbative corrections due to instantons, of order exp(— l//i). To understand these 
quantum-mechanical instantons in our problem, with a non-conventional Hamiltonian, we need 
a general, geometric approach to non-perturbative WKB expansions, like the one proposed in 
[51, 52]. In this approach, instanton contributions are obtained by looking at the complexified 
curve 



where H{q,p) is the Hamiltonian of the model. Perturbative WKB expansions are associated 
to periods of the above curve around "A- type" cycles, while non-perturbative corrections to the 
WKB method are associated to "B-type" cycles. In the case of ABJM theory, the complexified 
curve is identical to the spectral curve (4.44), after an appropriate choice of the variables. Its 

''After the first version of this paper was submitted, numerical and analytical studies of the function A{k) were 
performed in [50]. The ansatz (5.70) proposed here turns out to be incorrect. The results of [50] strongly suggest 
that the function A{k) can be obtained by resumming the so-called constant map contribution to the free energy, 
and expanding it around k — 0. The resulting series reproduces (5.68). We would like to thank the authors of 
[50] for informing us of their result prior to publication, which prompted us to correct a minor sign error in the 
calculation of (5.68). 



H(q,p)=E 



(5.72) 



-32- 



Riemann surface looks as shown in the Fig. 5. Let us introduce canonical coordinates Q, P 
related to the q, p coordinates as 



Q = q, 



p + q. 



(5.73) 



This preserves the symplectic form. The coordinate P is chosen so that it has no monodromy 
along the contour B. Then in the large E limit 



n(E) « ^ / PdQ ! 
The instanton contribution is of order 



1 4E^ 
\ol {{q,p) : IpI + |g| < 2E} 



exp 



h 



PdQ 



where in the large E limit 



^ PdQ = 2E ■ 47ri + ©(e"^^). 



(5.74) 



(5.75) 



(5.76) 



Here we used that, in the interior of the upper-right tube, P = p + q = 2E + 0{e~'^^) for some 
constant c, and that the monodromy of Q around the tube is 47ri. The above period can be 
computed exactly with the results of [8] since the behavior (5.76) fixes it completely: 



PdQ = -2ie^7r3F2 



1 1 1 



3 e 



2E 



= 8i7rS + 0(e-^^). 



2' 2' 2'"^' 2' 16 

2E\ 



i i 1 \ 

0'0'-2/ (5.77) 



In fact, after the identification (4.45), this period is equal to — 4^^, where Ag is the strong 
coupling instanton action computed in [8]. For large energy, (5.77) gives a contribution to the 
density of states of order 

exp [-4E/k] (5.78) 



which becomes a contribution 

exp [— 4/i/fc] 

to the grand canonical potential, and a contribution 



exp 



-2'K^2Njk 



(5.79) 



(5.80) 



to the canonical free energy. This is precisely the weight of a worldsheet instanton (2.6) in ABJM 
theory. 

Quantum-mechanical instantons are of course invisible in the perturbativc h expansion of 
and in the Wigner-Kirkwood expansion, but they appear in the 't Hooft expansion. In fact, 
the 't Hooft expansion of the canonical free energy 



F(A,fe) = ^fe2-2,^,(A) 

leads to a genus expansion of the grand canonical potential of the form [23] 

J't Hooft (/^, A;) = Y,}?-''3jg{iilk). 

g>0 



(5.81) 



(5.82) 
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Notice that, in the Fermi gas approach, only the perturbative part in n of J(/i) can be written in 
this form (5.82). The membrane instanton contributions and the function A(k) do not have the 
right functional dependence in /x/ A; to fit into the 't Hooft expansion, while the weight associated 
to a quantum-mechanical instanton (5.79) is again of the right form. In the case of ABJM theory, 
we see from (5.24) that the Fermi gas approach gives 

JgiO = O (e-^^) , g>2, 

where C = A*/^- From the point of view of the topological string, it follows from (2.30) that 
the variable C, is essentially the period T at large radius, and a perturbative Fermi gas approach 
makes possible to recover the leading, perturbative genus zero and genus one free energies of the 
topological string given in (2.25). 

Finally, we should mention that there is an extra source of worldsheet instanton-like correc- 
tions. In general, the exact representation (3.13) and the saddle-point integral (3.14) are only 
equivalent up to exponentially small corrections in N . Since we are taking into account such 
corrections, we have to be more careful here. The expression (3.13) is equivalent to 

Z{N) = — / d/x exp [J(/x) - /xiV] , (5.84) 

where the integration contour is parallel to the imaginary axis, and is arbitrary. To apply 
the saddle-point method, one chooses for /x* the saddle-point of the exponent, and then extends 
the integration contour to infinity along the imaginary axis (this is what gives the Airy function 
behavior we have found many times in this paper). As it is well-known, it is in this last step of 
extending the integration contour that one introduces exponentially small errors in A'^. A rough 
estimate of these errors can be done as follows. The saddle-point expansion involves integrating 
a Gaussian of the form 

-(/x-/x*)^j"(/x*) 



exp 



(5.85) 



The error in going to (3.14) can then be estimated by evaluating this Gaussian at the true 
endpoints in (5.84). This gives, by using the leading term in (5.24), 

~ exp[-2/i*/A;] (5.86) 

which is the square root of (5.79). Therefore, these type of corrections should also be taken in 
account when trying to extract information about worldsheet instantons. 



6. More general Chern— Simons— matter theories 

In this section we consider in detail more general CSM theories. We first study the thermody- 
namic limit of necklace quivers, and derive a general formula for the large N limit of their free 
energy which agrees with the result obtained in [17, 18] by analyzing the matrix model. Then we 
extend the considerations of section 5.3 to the general necklace CSM theories considered in sec- 
tion 3.2. For technical reasons we restrict ourselves to theories whose Hamiltonian is Hermitian, 
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i.e. whose free energy is real, and we show that, with that assumption, the Airy behavior of the 
resummed expansion found in [9] is indeed generic. These general considerations are then 
illustrated in detail in the case of the AB JM theory (i.e. the two-node theory) with fundamental 
matter. Finally, we consider the massive theories of [53], and we derive the N^^^ behavior found 
in [54, 14] with our techniques. 

6.1 Thermodynamic limit for general necklace quivers 

In this section we study the general necklace quiver considered in subsection 3.2. From (3.55) it 
follows that, for large energy, the Fermi surface is defined by the polygonal equation 



E 



(6.1) 



As an example, we show in Fig. 6 the classical Fermi surface (6.1) at small and large energy. 





-200 -100 



Figure 6: The Fermi surface for the three-node quiver with na = (1, 3, —4) in the q-p plane, for i? = 5 
(left) and E = 500 (right). At large energy it approaches the polygon (6.1). 

for a three-node quiver with = (1, 3, —4). By the by now familiar argument of the previous 
sections, the number of eigenstates is given by the semiclassical formula (4.13) applied to the 
domain bounded by (6.1). One finds. 



n{E) ^ CE"^, 



(6.2) 



where the constant C is given by 



Ki=l 



\x\ < 1 



(6.3) 



and the variables y,x differ from p,q in (6.1) by rescaling. Once this constant has been deter- 
mined, the large asymptotics of the free energy is given by (4.11) 



FiN) 



(6.4) 
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In order to compute C, we notice that the right hand side of (6.3) is the volume of a convex 
polygon which can be easily calculated. Suppose for simplicity that Nf. = 0. Let us introduce 
new parameters cj, related to ki in the following way 

Ca{j) =c + ^ki (6.5) 

i=l 

where c is an auxiliary constant and cr is a permutation chosen so that Cj < Cj+i,Vz. Then 



TT^C = Vol < 



{x,y) ■■ ^\y-Cjx\ <l \ ■ (6.6) 



Let us note that the Cj differ by a permutation from the parameters qa introduced in [10] (where 

they were defined in such a way that ka = Qa+i — Qa)- Notice also that the expression (6.6) is 
explicitly invariant under permutations of the cj. The inequality Yl^j=i Iv ~ < ^ defines a 
convex hull of 2r points (±Xs,=bys) so that 

r 

ys = CsXs, '^\ys - CjXs\ = I. (6.7) 
One finds ^ 

^ v^^^ — \ — z — I ' " v^^^ — f~z — I • ^^■^'> 

Then one can use the standard formula for the area of a convex hull to find 

TT^C = E '''^\7^'', V (6-9) 

where as usual we use the convention Cr+i = ci. Let us illustrate this formula by applying it 
to necklaces with three and four nodes. For the necklace with three nodes {k\,k2,kz) we can 
assume, without loss of generality, that 

|ci - C2I = l^al, |c2 - C3I = |ci - C3I = |A;2|- (6.10) 

Then 

= |fcl||fe2| + |fc2||fe3| + |fe3||fci| 

2 (|A:i| + |fc2|)(|fc2| + |A;3|)(|A;3| + |A;i|)- ^ ' ' 

For the quiver with four nodes, let us assume without loss of generality that Yli=i Ci = 0- Then 
an easy computation gives 

TT^C 1/1 1 ^ ^ 1 , ci + C3 



0.1 + 12 +4 , ' ' \, ]. (6.12) 

32 \Ci C4 C3 + C4 {cs + C4yj 

These formulae for the three and four-node quivers agree with the results first found in [10] 
(for the four- node quiver, their formula is obtained by setting ci = qs, C2 = qi, Cs = q2, C4 = q^). 
In fact, the above general result for the free energy of these CSM theories, involving the area 
of the polygon (6.3), has been derived in this form in [17, 18] by refining the analysis of the 
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matrix model done in [10] (where a different, but equivalent general formula for the free energy 
was proposed). This class of CSM theories is dual to M-theory on AdS4 x X^, where X-j is an 
appropriate tri-Sasaki Einstein space [39] . Therefore, the coefficient C should be proportional to 
the volume of the X'j manifold, and one should have 



Vol(X7) 
Vol(S7) 



— — = -\o\{ 
2 2 





x\<l 



> . 



(6.13) 



This is indeed the case, as it was proved in [17]. 

We then see that the Fermi gas approach allows us to rederive the result for the large A'' 
energy obtained in [17], but in a simpler way. The polygon appearing in the matrix model 
analysis of [17] has here a very simple interpretation: it is the large energy limit of the Fermi 
surface for the ideal Fermi gas. 

6.2 Airy function behavior for a class of CSM theories 

We now extend the considerations of section 5.3 to more general necklace quivers with matter. 
In the next subsection we apply the general considerations developed here to the case of ABJM 
theory with matter. 

Let us assume that the Wigncr transform of the density matrix can be written in the form 



where 



Qr = CLRq + bRp 



(6.14) 



(6.15) 



for suitable a^, b^, and the different Qr are given by linearly independent combinations of q,p. 
We will also suppose that the functions $j are real valued, even, and that^ 



$i(Q) = 7i|Q|+C?(e-^l«l), Q^oo, 7i>0. 



(6.16) 



These assumptions are obviously true for the general necklace quivers considered in subsection 
3.2. In addition, we will suppose that H\\i{q,p) is real. This corresponds to the case when the 
quantum operator H is Hermitian. It should be possible to treat the general case, when H has 
complex eigenvalues, with similar techniques, but we will not do it here. 

As usual, the leading contribution to the number of states is given by the volume of a polygon: 



niE) « ^Vol |(g,p) : f^-filaR^q + bR^p\ < = CE\ 



(6.17) 



Let 



CE = {{q,p) : H^{p,q)=E} (6.18) 

be the real curve describing the Fermi surface. One can always decompose Ce into patches Ur, 
U'ji so that both Ur and U'^ contain one point where Qr = 0, Wjj is related by to Ur by reflection 
through the center of the polygon, and 



Yli\aRiq + bRiP\ 



E + 0{e 



-cE\ 



(6.19) 



dUR,au' 



*In what follows, the constant c > may have difTerent values in different formulae. 
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Figure 7: The black thick Une depicts the Fermi surface (6.18). The green thin line depicts the limiting 
polygon. The red thick dashes mark the boundaries of the patches Un. 



In Fig. 7 we depict the general structure of such decomposition. In the case of ABJM theory 
considered in section 5.3, the boundaries of the regions I and II lying on the curve Ce are halves 
of the patches lAi and hl2- A particular example of such a decomposition of the Fermi surface 
in regions is shown in Fig. 10, in the case of ABJM theory with fundamental matter. We will 
now argue that, for each patch Ur, the structure of corrections is essentially the same as in the 
ABJM case. 

One can always choose 

Pr = CRP + dRQ (6.20) 

so that 

dQR A dPR = dg A dp (6.21) 
and such that, in the domain Ur, the Hamiltonian can be written as follows: 

H^{q,p)=/3RPn+ '^^{QR)+Y1 E ^'''c^M''\QR)+0{e-'^'') (6.22) 

i\Ri=R r>Oi\R^=R 

where 0{e~'^^) denotes an estimate which is uniform in Ur. Without loss of generality one can 
assume that Pr > 0. Then in the domain Ur the solution to H^y [Pji, Qji) = Hw{p, l) = E can 
be written as 

Pr{E, QR) = j-{E- r^\QR\ I + ^Pr{E, Qr). (6.23) 

V i\R,=R ) 

In this equation, 

/^Pr{E, Qr) = ApPniQn) + 0{e-^^), (6.24) 
where Ap denotes the perturbative part of the correction. It is given by 

A^Pr{Qr) = I E (^^(^«) - + E E ^'''c^M''\Qr) ] ' (6-25) 

\i\Ri=R r>Oi\R^=R J 
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and satisfies the property, 

^pPr{Qr) = O (e-^l'^^l) ,Qr^oo (6.26) 
The property (6.19) imphes that 



Vol{Hw{q,p) <E} = Vol <j Y.^i\QR\ < E \ + 2 ^ AVoli^(£;) + C(e-'=^) (6.27) 

where 



, i=i ) R 



From (6.24) it follows that 



AVoIr{E) = J APR{E,QR)dQR. (6.28) 



AYolRiE) = J A^Pr{Qr)AQr + C»(e-^^). (6.29) 



By using (6.26), we can extend the integration region to infinity, up to non-perturbative correc- 
tions, and we obtain 



AYo\r{E)= j ApPR(gR)dQR + 0(e-'=^). (6.30) 



Let us denote 



oo oo 

\Pr{Qr)AQr = -^ / {'^^iQR)-l^\QRWR-^ E (6-31) 

PR J ., J-, PR ., 

— oo — oo — ^ — -rt 

Similarly to what happened in ABJM theory, there are no perturbative corrections to n{E) from 
higher terms of the Wigner-Kirkwood expansion. Therefore we obtain 

n{E) = CE^ + no + Oie'"^) (6.32) 

with 

rio = ^^ApVolfl. (6.33) 

R 

As was shown in subsection 5.2, it then follows that the 1/N corrections are resummed to an 
Airy function 

Z{N) = C-^^^e^ Ai \c-^/^{N - B)] + Z^piN) (6.34) 

where 

5 = no + — . (6.35) 

The above general argument gives an explicit algorithm to compute the constant (6.33). In the 
next subsection we consider the example of the ABJM theory with matter as an illustration of 
this argument. 
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Figure 8: Quiver for the two-node theory with fundamental matter. 



6.3 ABJM theory with fundamental matter 

The ABJM theory with matter we will consider is described by a two-node quiver with equal 
number of fundamentals in each node, see Fig. 8. The density matrix of this theory is given by 
(3.57): 

pw = exp J-i7w} = e-^(«)/' * e'^^"^ * e'^^^ * e'^^"^ * e-^(«)/2 (6.36) 
where U{q), T{p) are given respectively in (3.30) and (3.31), and 

p + q 



*(P + q) = Nf log 2 cosh - 



2k 



(6.37) 



The Wigner transform of the Hamiltonian has the following H expansion, which can be obtained 
with the use of the Baker-Campbell-Hausdorff formula: 

H^ip, q) = Ha{p, q) + Y^U'{qf {2^"{p + q) + T"{p)) 



12 



(-2V(p + qf {T"{p) - 2U"{q)) + T'{pf (2V'(p + q) + U"{q)) + 4T' {p)U" {q)V {p + q)) 



+ -U'{q) {T"ip)^'{p + q)- T'ip)^"{p + q)) + O {h^) , 

where the first term is the "classical" Hamiltonian^: 

Hd{p, q) = T{p) + U{q) + 2^{p + q). 

The function ^ has an asymptotic behavior similar to that of T and U: 

a , 



where 



N 



a 



f 



k 



(6.38) 
(6.39) 

(6.40) 
(6.41) 



Since i^w is real, (6.36) is a particular example of the general case (6.14) considered in the 
previous subsection. We have three different local coordinates 



Qi = q, Q2=P, Q3=P + q- 
For large energy the Fermi surface H^^{q,p) = E approaches a polygon given by 

^4+^4+a\p + q\=E, 



(6.42) 



(6.43) 



^We use quotation marks because * still contains k = h/{2Tv). 
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Figure 9: The Fermi surface (6.18) in the q-p plane for the ABJM theory with fundamental matter and 
TV/ = 1, fc = 2, for £; = 5 (left) and E = 100 (right). 



see Fig. 9. Therefore, the leading contribution to the number of states is 

""^^^ 2^^°^ I ^ + ^ + a|p + gl < = Ci?2 (6.44) 

where 

C= + (6.45) 

It follows that 

3 Vl + a 

which reproduces the result of [11, 10]. Notice that, as in the case of ABJM theory, the large 
energy limit of the Fermi surface is closely related to the tropical limit of the spectral curve 
obtained in [11]. 

Now let us compute the corrections according to the general scheme described in the previous 
subsection. The regions UrM'r for i? = 1, 2, 3, as well as the lines Qr = 0, are shown in Fig. 10. 
In the domain Ui the Hamiltonian can be written as 

^w(5, P) = P (a + ^) + U{q) +qa- ^h\2a + lfU"{q) + ^i^n^i ^(2n) ^ ^(^-cS)^ (g ^^) 

^ ^ n>l 

Therefore we can take 

Oi 

^i=^+T7^^« (6.48) 
1/2 + a 

and 

I n>l ) 
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Figure 10: The regions Ur, U'^ defined in section 6.2 for the Fermi surface of ABJM theory with matter. 
The dashed hnes are defined by Qn. = 0, where the different coordinates are given in (6.42). 



In the domain U2 we have, 
Hwiq, P) = Tip) -pa-q(^a + ^^+^i2a + lfh^T"{p) + ^ h'^clT^^^^ (p) + ©(e"^^), (6.51) 

therefore 



n>l 



and 

ApP2(Q2) 



P2 = -q 



a 



1/2 + a 



p 



(6.52) 



l + 2a 



+ ^^h\2a + lfT"{Q2) + ^'"c2t(2-)(Q2) 1 , (6.53) 

n>l J 



ApVob 



3(1 + 2a) 48 



(6.54) 



In the domain , 

^w(g,p) = p/2 - q/2 + 2^>{p + q) + ^hH" {p + q) + Y ^i^n^s ^(2n) ^ ^) ^ ©(e"^^). (6.55) 

48 



n>l 



Therefore 



and 



p-q 



(6.56) 



ApP3(Q3) = - {2^m -am)- ^hH" (Q3) - E ?i'"4xI/(2n)(Q3), (6.57) 

(6.58) 



n>l 



ApVols = -^^^a 
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Finally, we obtain 



1 



J2 ApVolij 



(6.59) 



•Re{i,2,3} 

and we conclude that the partition function is given by the Airy function 

Z{N) = C-^/^e^ Ai \c-^/^{N - B)] + Znp(iV) 

with 



B = no + 



Ctt' 



1 



+ (1 - 3a) . 



(6.60) 



(6.61) 



3 3A;(1 + 2a)2 24 

The above procedure can be repeated for other quivers with a Hermitian Hamiltonian in 
order to determine the precise value of B. For example, for the four-node quiver with levels 

{k,-2k,2k,-k), (6.62) 

one finds 



B = - 



13 



+ 



k 



135k 8 



(6.63) 



Clearly, it would be nice to have a close answer for the shift for a more general class of quivers 
(like for example four node quivers with a Hermitian density matrix). In addition, it would be 
interesting to compare the shifts (6.61), (6.63) with a direct calculation from the M-theory/type 
IIA geometry, as in [27, 28]. 

6.4 The massive theory 

The techniques developed in this paper can be also applied to a variant of ABJM theory in which 
the Chern-Simons levels ki, k2 do not add up to zero [53]. We will denote 



.^11 47r 1 

27ri^ = ----, -Tr = Tr 
ki k2 n k\ 



1 

k2' 



(6.64) 



so that the original ABJM theory is recovered when ^ = 0. Notice that 6 is in principle imaginary, 
but it will be useful to Wick-rotate it to real values (sec also [56]). The theory with ki + k2 ^ ^ 
was studied in [53], where it was argued that a non-zero 9 corresponds to a non-zero Romans 
mass in type IIA supergravity. For this reason, we will call this theory the "massive" theory. 
The massive theory was further investigated in [54] , where it was found that its free energy scales 
with TV as 

F{N) « (fci + k2f'^ N^'^. (6.65) 

This scaling was reproduced in [14] from an analysis of the matrix model representing the partition 
function, 



z{N,e) 



d^x d^^ y 



n 



2 sinh 



2 



2 sinh 



(27r)^ (27r) 



N 



n 



2 cosh 



exp 



N 



(6.66) 



1=1 
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The exact planar resolvent of this theory was found, in a somewhat implicit form, in [55]. The 
scaling (6.65) can be also derived from this resolvent by using the techniques of [56]. 

In order to apply the Fermi gas picture to this theory, we have to find an appropriate density 
matrix. An elementary computation leads to 



aeSN 



where 



p(xi,X2;9) = e-^t^«(^i)i^(xi,x2;e)e-^^«(^2)_ 
Here, the one-body potential is given by 



Ue {q) = log ^2 cosh 



while the function K is given by 
K{xi,X2;0) 



dy 



^ AttH cosh I 



exp 



[V^ - y 



-{xi +X2) + -^{xi - X2) 



Although (6.70) is complicated, its Wigner transform is very simple, 
where 

P\ 



Tg{q,p) = log (2 cosh + -p2 + Qpq^ 



and of course 



The Wigner transform of the density matrix is then 

pw(^) = ^/lTWWJA e-i^''('?) ★e-^^^^'^') ^e-i^''^"), 
and defines the Hamiltonian of the theory through 

pw(^)=e.-''*('^ 



(6.67) 



(6.68) 



(6.69) 



(6.70) 



(6.71) 

(6.72) 
(6.73) 

(6.74) 
(6.75) 



For 9 = we recover the density matrix of ABJM theory (3.39). 

We can now use the technology developed before to analyze the theory. We will content 
ourselves with an analysis of the thermodynamic limit, which leads to a nice interpretation of 
the N^/^ behavior found in [54, 14]. We will also assume that 



or equivalently, that 



\eh\ < 1 

ki + k2 



fei - ko 



< 1. 



(6.76) 
(6.77) 
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Figure 11: The Fermi surface (6.81) in tlie q-P plane, for = 1/10, k = 1, E ^ 50 (left) and E = 5000 
(right). When the energy is large, the Fermi surface approaches the surface defined by (6.82). 



In this limit we can safely ignore quantum corrections and just look at the classical Hamiltonian 
H,i{q, P; 9) = U{q) + T {P - q) + ^P^ - ^ log (l + ^) , (6.78) 



where 

P = p + q. (6.79) 
This linear change of variables preserves the volume form in phase space, 

dg A dP = dg A dp. (6.80) 

At large E the Fermi surface 

H,iiq,P;9)=E (6.81) 

becomes simply 

%2 + = E, (6.82) 



as we can see in Fig. 11. Notice that, once 6^0, the equation defining the Fermi surface at 
large E has a quadratic term in the new momentum coordinate P which dominates at large E. 
In other words, the Fermi gas has now a non-relativistic dispersion relation, and this changes the 
scaling of the free energy. Looking at (4.15) we deduce that 

s = l, (6.83) 
therefore the free energy should scale now as iV^/^ as found in [54, 14]^. We find 

/ 2E I — 

n(E) ^A^r'dPiE- ^-pA = -^JIe^'. (6.84) 
^ ' 2-Kh L I 2 / 37r/i V (9 ^ ^ 



^The matrix model analyzed in [23] displays tlie same scaling. 
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The free energy can now be computed from (4.11) and reads, 



F^N) ^ (^^^y%^/^N'/\ (6.85) 
If we express this in terms of ki + k2 we find, 

q5/3 / /i2fc2\l/3 

Since we are assuming (6.76), our result can be written as 

o5/3 

F{N) ^ -—^ne-f{k, + k2)'/^N^/^ (6.87) 

which is precisely what [14] obtained. Notice that, in [14], this resuh was derived based on an 
assumption on the behavior of the eigenvalues of the matrix model at large N, while here we 
have obtained it directly. When the parameter d^Ji^ is not small, one has to take into account 
the quantum corrections to the Hamiltonian, and the equation of the Fermi surface is modified. 
It would be interesting to study in more detail the different regimes that can occur in this theory 
as we vary the coupling constants. 



7. Conclusions and prospects for future work 

In this paper we have developed an ideal Fermi gas approach to a large class of CSM theories with 
> 3 supersymmctry. Since the particles in the gas do not interact, all the information of the 
problem is encoded in the one-particle quantum Hamiltonian. We have seen that the structure of 
the CSM theory determines the detailed form of this Hamiltonian, which is conveniently encoded 
in its Wigner transform Ht^, 

CSM theory — > quantum Hamiltonian -ffw- (7-1) 

The semiclassical analysis of this Fermi gas already gives us a lot of information on the original 
partition function. As we have seen, the leading semiclassical contribution provides an extremely 
simple derivation of the N^^"^ behavior in these theories (including the right coefficient), and the 
next-to- leading correction already encodes the full 1/N expansion of the original matrix model. 
The large energy limit of the Fermi surface is a polytope which encodes the geometry of the 
dual tri-Sasakian target geometry. In addition, non-perturbative membrane instantons can be 
computed order by order in k (but at all orders in the membrane winding). We summarize some 
of the features of the Fermi gas approach, as compared to the 't Hooft expansion, in the following 
table: 





't Hooft expansion 


Fermi gas 


semiclassical limit 


spectral curve 


Fermi surface 


quantum corrections 


quantum spectral curve 


quantum Hamiltonian 


exponentially small corrections 


worldsheet instantons 


membrane instantons 


non-perturbative effects 


membrane instantons 


worldsheet instantons 
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It would be interesting to understand at a deeper level the Fermi gas picture. Since D-branes 
behave as fermions (see for example [57]), and the gauge theories we have considered have D- 
brane realizations, one might be able to derive this picture directly from the D-branes underlying 
the gauge theory. 

Our formalism has many similarities with recent developments matrix models and topological 
string theory (see, for some examples, [57, 58, 59]). The semiclassical limit gives a Fermi surface, 
which can be regarded as the counterpart to the spectral curve in the conventional 't Hooft 
expansion. For ABJM theory both curves are the same, but for the more general CSM theories 
studied in this paper, the semiclassical Fermi surface is simpler and easier to find. In fact, the 
planar resolvent of the CSM matrix models is in general not known, and when it is known (as in 
[55, 11]), it turns out to be a complicated transcendental function. It seems however that both, 
the Fermi surface and the spectral curve, have the same tropical or polygonal limit. 

The Fermi surface has quantum corrections, and one then obtains a "quantum corrected 
curve" nicely encoded in the function ifw in phase space. The quantum corrections appear in a 
very natural way, by replacing standard products by * products. However, the corrections found 
in this way are not in the string coupling constant gs, as in conventional topological string theory 
[57, 60] but rather in the inverse string coupling 1/gs- Therefore, the quantum Hamiltonian 
is computing the quantum curve for the strongly coupled string. 

More generally, our formalism seems appropriate to study the strong-coupling regime of 
matrix models and topological strings. For example, in the case of ABJM theory, it is known 
[26, 5] that the ABJM matrix model (2.1) corresponds to the submanifold Ti = T2 of the moduli 
space of topological strings on local x P^. Here, Ti^2 are the Kahlcr parameters measuring the 
sizes of the two P^s in the geometry. However, by looking at for example (5.83), it is clear that the 
grand canonical potential of ABJM theory seems to be directly related to the topological string 
free energy in the large radius frame. Therefore, our calculation of J(/u) for the ABJM theory 
can be interpreted as a concrete strong coupling expansion of this topological string free energy, 
including non-perturbative effects. The worldsheet instantons of the topological string at large 
radius would appear then as quantum-mechanical instantons of the Fermi gas. Notice also that 
the grand canonical partition function, which is the focus of this paper, involves the sum over 
fluxes first considered in the context of topological strings in [57], and studied from the matrix 
model point of view in [61, 62, 63]. Our formalism gives a concrete approach to calculate this 
object at strong coupling, but one should clarify the relation between the picture proposed here 
and the non-perturbative approach of [61, 62] involving theta functions on the spectral curve. It 
would be very interesting to develop further all these relationships to topological string theory. 

If we stay in the world of CSM theories and their AdS duals, an obvious and important 
question is to which extent the Fermi gas formalism can be applied to other theories. For 
example, one should consider, already in the = 3 case, the issue of quantum corrections and 
possible Airy behavior in theories where the free energy is not real, as well as in theories where 
the nodes have different ranks [64]. The extension to A/" = 2 theories is more challenging. In that 
case the interaction among eigenvalues involves a more complicated function [3, 4] and they lead 
generically to interacting Fermi gases, rather than to ideal gases. 

Although the equivalence between the matrix models and the Fermi gas partition function is 
an exact statement at finite N and k, in this paper we have worked in the thermodynamic limit 
(large N) and in the semiclassical limit (expansion in powers of k around A; = 0). Fortunately, 
as we have seen, the expansion in l/N satisfies some sort of "non-renormalization" property and 
we can determine it at finite A; by a next-to-leading computation in the WKB expansion. An 
obvious challenge is to solve the Fermi gas problem at finite k, say A = 1, to make full contact 
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with the M-theory expansion. This would amount to a resummation of the non-perturbative 
effects computed in this paper order by order in k. One possible route to achieve this is to find 
the exact eigenvalue spectrum of the quantum Hamiltonian, or equivalently, of the density matrix 
p. In the case of ABJM theory, this means solving the integral equation (3.19), which can be 
written as 

J dx'p{x,x')(j)n{x') = XnMx), (7-2) 
where the kernel can be written in the form 

p{x,x') = , . (7.3) 

47rA;cosh (^2^) 

This type of kernel appears in other contexts, like the 0(2) matrix model [22] (albeit with a dif- 
ferent U{x)), and it is connected to both the Hirota hierarchy [22, 65] and to the Thermodynamic 
Bethe Ansatz [66, 67]. At least for = 1,2 (where supersymmetry is enhanced), we anticipate 
a nice solution to the eigenvalue problem in terms of an integrable system. In particular, the 
relation to differentiable hierarchies of the Hirota type suggests that the conjecture 1.1 might be 
proved by performing a suitable double-scaling limit in the hierarchy. 

It would be also very interesting to develop further the relationship between worldsheet 
instantons and quantum mechanical instantons of the Fermi gas sketched in section 5.4. It is 
clear that the solution of the ABJM theory found in [5] , in the 't Hooft expansion, is extremely 
powerful in order to capture these corrections, but it would be important for the development of 
the Fermi gas approach to have a better understanding of this issue. 

Finally, there might be a connection between the Fermi gas of this paper and two other 
pictures for string/M-theory based on fermions: the droplet picture proposed in [68, 69] to 
analyze 1/2 BPS operators in A/" = 4 SYM theory, and the Fermi liquid picture of non-critical 
M-theory proposed in [70]. 
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A. Quantum corrections in ABJM theory at order 0{h'^) 

In this Appendix, we give some details on the computation of the order 0{h'^) corrections to the 
grand canonical potential of ABJM theory, which confirm the general arguments of section 5. 
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The Baker-Campbell-Hausdorff formula applied to (3.39) gives 
H^{q,p) = T + U + ^[T,[T, UH + ^ [U, [T, UH + ^[[[[T, [/]., U]., U],, T], 

- ^ im, TU U],, T],, U], + ^ [[[[[/, r]„ T]„ T],, U], + ^ [[[[T, T],, T], 

+ 5^ ^ttf^' " 7^ ftft^' + ■ ■ ■ 

This leads to the next correction to the Wigner transform of the Hamiltonian 

1 ^\2rW// V /„/// s\2 , 1 /ml/ \\2 mil / \ l-nll/ \\2 



(A.l) 



1 (c/'(g))V"(g)T'(p)r'"(p) + ^ (r'(i,))V"(p)c/'(g)c/'"(g) 



(A.2) 



The computation of J2(/u) also involves the defined in (5.3) up to order 0(%^\ Due to (5.5) 
only the terms with r < 6 are needed. A long but straightforward calculation leads finally to 

^^^^^ = "4320 " 2880 ^^"^ + ^""^ " ^'^^^ + ^^^^^ '""^ + ^ ^^^""''^^ ' ^^'^^ 
Notice that no polynomial in is generated, as expected from the analysis in section 5. 
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